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PREFACE 

This  is  the  second  in  a series  of  reports  on  sene  descriptive  studies 
of  cyclones  which  were  carried  out  as  part  of  a research  project  on  cyclo- 
genesis (Project  SCUD)  sponsored  by  the  Office  of  Naval  Research.  (The 
first  paper  in  the  series,  by  E.  L.  Fisher,  dealt  with  "Upper  level  temper- 
ature fields  associated  with  east  coastal  cyclogenesis".) 

There  exists  a need,  both  by  the  dynami cists,  who  are  attempting  the 
monumental  task  of  weather  forecasting  by  hydrodynaroical  methods,  and  the 
empirical  forecasters,  for  quantitative  information  about  the  structure  of 
cyclone  development.  Whether  such  information  is  used  to  construct  new 
synoptic  models,  to  provide  parameters  for  statistical  prediction  schemes, 
or  to  suggest  methods  for  solving  the  development  problem  in  numerical 
forecasting,  its  essential  importance  is  obvious. 

One  device  for  summarizing  synoptically  many  facts  about  the  structure 
of  some  atmospheric  phenomenon  is  the  composite  map.  Composite  maps  have 
been  used  rather  rarely  in  meteorology.  let  they  do  provide  a convenient 
way  to  test  old  models  of  the  atmosphere  and  to  construct  new  ones.  In  the 
first  section  of  this  report  the  composite  map  technique  has  been  employed 
to  construct  some  pictures  of  the  average  horizontal  temperature  distribu- 
tion in  the  lower  stratosphere  over  cyclones  which  do  and  cyclones  which  do 
not  produce  secondary  lows  along  the  east  coast  of  the  United  States.  This 
study  in  comparative  analysis  was  undertaken  by  Mr.  Emanuel  K.  Ballenzweig, 
Research  Assistant,  to  test  the  hypothesis  of  Mr.  E.  L.  Fisher  that  there 
may  be  found  a characteristic  anomaly  in  the  200  mb  temperature  field  prior 
to  the  development  of  certain  types  of  east  coastal  secondary  cyclones.  (The 
article  by  Mr.  Ballenzweig  is  part  of  an  M.S.  thesis  submitted  to  the  Depart- 
ment of  Meteorology  and  Oceanography  at  New  York  University.) 

The  routine  analyses  of  upper  level  maps  reveal  that  deep  upper  level 
cyclones,  often  with  no  corresponding  sea  level  cyclone,  occur  frequently 
and  play  an  important  role  in  low  level  cyclogenesis.  These  systems  are 
characterized  by  cold  tropospheric  cores  and  low  warm  stratospheric  cores. 

The  development  and  propagation  of  the  thermal  fields  in  these  cyclones 
must  be  affected  to  an  important  degree  by  vertical  motions  in  the  atmosphere. 
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The  second  article  in  this  report  presents  some  calculations  of  the  vertical 
velocity  fields  at  500  and  200  mb  in  such  a deep  cold  cyclone.  This  research 
was  done  by  Lieutenants  Benjamin  M.  Heimuu  and  John  J.  R.  Kinney,  U.  S.  Air 
Force,  while  on  temporary  duty  at  New  York  University.  (Although  they  were 
not  employed  on  the  project,  Lieutenants  Herman  and  Kinney  worked  in  close 
collaboration  with  Project  SCUD  and  have  authorised  the  inclusion  of  an 
abridged  version  of  their  k.S.  thesis  in  this  report.) 

Interest  in  the  analysis  of  temperature  fields  at  upper  levels  has  been 
stimulated  recently  both  by  the  revival  of  certain  thermal  theories  of  cyclo- 
genesis as  well  as  by  theoretical  developments  in  the  use  of  baroclinic 
models  in  numerical  forecasting.  For  a long  time  it  has  been  customary  in 
synoptic  meteorology  to  assume  that  radiative  temperature  changes  are  neglig- 
ible except  within  a few  thousand  feet  off  the  ground  and  perhaps  in  the 
ozonosphere.  Such  an  assumption  (especially  when  coupled  with  the  assump- 
tion that  temperature  is  advected  horizontally)  leads  to  great  simplification 
in  theory  and  practice.  However,  it  is  of  some  importance  to  determine  the 
magnitude  of  this  neglected  quantity.  Some  calculations  of  radiative  temper- 
ature changes  in  cyclones  are  presented  in  the  paper  by  Mr.  George  Ohring, 
Research  Assistant.  These  calculations  are  helpful  for  the  understanding 
of  certain  transient  characteristics  of  cyclonic  temperature  fields,  espec- 
ially near  cloud  tops,  and  in  the  interpretation  of  lapse  rates.  (The  paper 
by  Mr.  Ohriug  is  an  abridgment  of  an  M.S.  thesis  presented  to  the  Department 
of  Meteorology  and  Oceanography  at  New  York  University.) 

The  last  article  in  this  report  contains  a supplement  to  James  E.  Miller's 
statistical  study  of  east  coastal  cyclogenesis. 


Jerome  Spar 
Project  Director 
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COMPOSITE  200  MB  TEMPERA TUH£  FIELDS  ABOVE 
E1TRATR0PICAL  CYCLONES 


EMANUEL  M.  fa ALLEN ZWEI G 

Introduction 

The  suggestion  has  often  been  made  that  cyclogenesis  is  the  result  of 
high  level  density  changes.  Although  this  view  is  apparently  supported  by 
the  well  known  observation  that  deep  extra tropical  cyclones  are  character- 
ized ty  warm  stratospheres,  the  causal  relationship  between  the  cyclone  and 

the  warm  stratosphere  is  still  uncertain, 
n 

Stuve  (1927)  proposed  toe  theory  that  cyclogenesis  and  anti cyclogene sis 
are  the  result  of  meridional  air  motions  in  the  stratosphere.  In  toe  tropo- 
sphere too  temperature  gradient  iB  negative  from  the  equator  to  toe  poles; 

II 

in  the  stratosphere  the  gradient  is  positive.  According  to  Stuve,  a surface 
pressure  minimum  is  produced  ly  the  equatorward  transport  of  warm,  light 
stratospheric  air  over  warm,  light  tropcBpheric  air.  This  model  stresses 
the  stratospheric  density  changes,  but  neglects  toe  compensating  tropo- 
spheric density  changes  which  were  observed  ty  Dines  (1919). 

Wulf  and  Obloy  (1944)  have  also  hypothesized  that  surface  pressure 
changes  are  due  to  the  interaction  of  stratospheric  end  tropospheric  temper- 
ature changes.  They  visualize  a warm  artic  stratospnere  undercut  by  a warm 
tropical  troposphere.  According  to  their  toeoiy  , this  combination  produces 
a pressure  minimum;  frit  this  theory  also  does  not  quite  account  for  the  com- 
pensating effect  of  the  tropospheric  density  changes. 

The  thermal  theories  of  cyclogenesis  (e.g.,  Austin  1952)  are  based  large- 
ly on  the  principle  of  horizontal  advection.  However,  it  is  well  known  from 
studies  by  Fleagle  (1947)  and  others  that,  especially  in  the  stratosphere, 
vertical  motions  play  an  important  part  in  toe  development  of  temperature 
fields.  The  role  of  vertical  motion  is  clearly  illustrated  in  a recent  naper 
by  Fisher  (1954)  Bhowing  virtually  stationary  temperature  fields  at  200  mb 
embedded  in  rapidly  moving  air  currents. 
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Although  the  exact  causal  connection  between  surface  pressure  changes  and 
stratospheric  temperature  changes  is  not  known,  there  is  reason  to  believe  that 
the  stratospheric  temperaoure  fields  may  provide  indications  of  surface  devel- 
opments. Fisher  (1954)  has  suggested  that  there  may  be  a characteristic  anomaly 
in  the  form  of  a tongue  of  warm  air  in  the  200  mb  temperature  field  prior  to 
the  development  of  a secondary  cyclone  along  the  east  coast  of  the  United  States. 

This  investigation  was  designed  to  study  the  200  mb  temperature  field  above 
extratropical  cyclones  in  the  United  States  and  to  determine 

(1)  whether  there  is  any  difference  between  the  200  mb  temperature 
fields  and  temperature  change  fields  of  those  cyclones  which 
develop  secondaries  and  those  which  do  not, 

(2)  the  average  structure  of  the  200  mb  temperature  field  above 
moving  extratropical  cyclones,  and 

(3)  the  200  mb  temperature  field  above  the  newly  formed  cyclone. 


The  Composite  Method 

A study  was  made  of  ~ix  cyclones  which  formed  secondaries  to  the  south- 
east of  the  primary  cyclone  and  of  seven  cyclones  which  did  not  form  secondar- 
ies during  the  period  December  1952  - February  1954.  (see  table  1).  Those 
cyclones  which  formed  secondaries  are  called  group  1 cyclones;  those  cyclones 
which  did  not  form  secondaries  are  called  group  2 cyclones. 

The  cyclones  were  tracked  from  the  time  they  entered,  or  formed  within, 
o o o 

the  area  bounded  by  30  N latitude,  50  N latitude,  105  ¥f  longitude  and  the  east 
coast  of  the  United  States.  (It  was  considered  inadvisable  to  study  cyclones 
in  the  vicinity  of  the  Rocky  Mountains  because  large  errors  may  result  from 
the  UBe  of  nonrepresentative  surface  temperatures  in  making  pressure  reductions 
to  sea  level.  This  difficulty  is  almost  entirely  eliminated  east  of  the  Rockies.) 
Three-hourly  positions  of  the  cyclones  were  obtained  from  the  United  States 
Weather  Bureau  maps  at  LaGuardia  Field  after  they  were  detected  on  maps  at 
New  York  University.  These  positions  and  intensities  were  transferred  to 
blank  maps.  Tracks  were  drawn  using  only  the  0330Z  and  the  1530Z  positions. 

A standard  reference  time  (zero  time)  was  needed  in  order  to  compare  the 
cyclones  which  formed  secondaries  with  those  which  did  not.  _cro  time  for 
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those  cyclones  which  formed  secondaries  was  easily  defined  as  the  time  closest 
to  the  RAOB  time  at  which  the  secondary  formed.  However,  there  is  do  corres- 
ponding zero  time  for  those  cyclones  which  do  not  produce  secondaries.  The 
zero  time  for  the  latter  group  was  defined  in  the  following  manner. 

A map  (fig.  1)  was  constructed  showing  the  frequency  distribution  of  the 
locations  of  primary  cyclones  at  the  time  secondary  cyclones  first  appeared. 

This  map  was  based  on  the  data  collected  by  Miller  (1945)  for  type  B or 
warm  frontal  secondaries  which  form  to  the  southeast  of  the  primary  cyclone 
in  the  eastern  United  States.  This  is  the  type  of  development  with  which  the 
present  study  is  concerned = 

The  tracks  of  the  group  2 primary  cyclones  were  superimposed  on  figure  1. 
The  position  along  the  track  at  which  the  highest  frequency  value  was  found 
was  designated  as  the  zero  position,  and  the  corresponding  time  was  consider- 
ed to  be  the  zero  time.  The  zero  time  for  this  group  of  cyclones  is  thus  the 
time  at  which  a secondary  is  most  likely  to  appear  (although  none  did  develop). 

The  method  outlined  above  was  tested  on  the  six  cyclones  in  group  1 which 
did  produce  secondaries.  In  five  of  the  six  cases  the  zero  time  obtained  from 
figure  1 corresponded  to  the  observed  zero  time. 

« 

The  zero  times  for  the  cyclone*  studied  are  shown  in  table  1. 


Table  1.  Zero  time  of  primary  cyclones. 


Group  1 

Group  2 

03  Z 

21  December  1952 

15  Z 

24  January  1953 

03  Z 

11  February  1953 

15  Z 

16  February  1953 

15  Z 

12  February  1953 

03  Z 

21  February  1953 

15  Z 

23  March  1953 

03  Z 

4 March  1953 

03  Z 

11  April  1953 

03  Z 

8 March  1953 

15  Z 

28  October  1953 

15  Z 

4 December  1953 

15  Z 

26  February  1954 

In  comparing  the  two  groups  of  cyclones  from  the  time  they  entered,  or 
formed,  within  the  area  defined  above,  it  was  not  possible  to  go  back  in  time 
more  than  twenty-four  hours  from  zero  hour  nor  ahead  more  than  twelve  hours 
because  the  different  cyclones  had  different  life  spans  in  the  area  as  can  be 
seen  in  table  2. 


Fig.  1.  Isopleths  of  frequency  of  location  of  primary  cyclones 

t 

at  the  time  of  secondary  cyclogenesis.  j 

! 
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Even  the  composite  maps  for  group  1 at  zero  plus  twelve  hours  and  for 
group  2 at  zero  minus  twenty-four  hours  are  of  dubious  validity  since  they 
contain  only  about  one-half  the  total  number  of  cases. 

Table  2.  The  number  of  primary  cyclones  in  each 
group  at  each  time. 


Group 

0-48 

0-36 

0-24 

0-12 

0 

0+12 

0+24 

hours 

1 

3 

4 

6 

6 

6 

3 

0 

2 

0 

1 

4 

7 

7 

7 

3 

The  following  procedure  was  used  to  construct  the  composite  maps: 

1)  The  surface  positions  and  central  pressures  of  the  cyclones  were  en- 
tered on  the  corresponding  200  mb  charts. 

2)  A grid  of  2°  squares,  13°  in  latitude  and  2 4°  in  longitude,  was  placed 
on  the  200  mb  chart  in  such  a way  that  the  horizontal  lines  of  the  grid  were  on 
the  odd  latitude  lines  of  the  map,  the  vertical  lines  of  the  grid  were  on  the 
even  longitude  lines  of  the  map,  and  the  cyclone  was  in  the  square  marked  X 
(see  figure  2) . 

3)  With  the  cyclone  center  as  the  ref ax  fence  point  in  each  case,  the 
temperatures  at  the  intersections  of  the  horizontal  and  vertical  lines  of 
the  grid  were  interpolated  to  the  nearest  0.5°C  and  plotted  on  a corres- 
ponding grid. 

4)  The  mean  temperature  at  each  grid  point  was  determined  and  a set  of 
200  mb  mean  temperature  charts  was  constructed  for  each  group  of  cyclones.  (A 
set  consists  of  four  charts:  zero  - 24  hours,  zero  - 12  hours,  zero  hour  and 
zero  + 12  hours) . 

f»)  Two  cots  of  200  mb  mean  temperature  change  maps  were  also  constructed, 
one  for  each  group  of  cyclones.  (A  set  consists  of  three  charts:  zero-24  to 
zero-12  hours,  zero-12  hours  to  zero  hour  and  zero  hour  to  zero  •+  12  hours) . 

6)  The  grid  was  also  used  to  study  the  200  mb  temperature  disti’ibution 
above  the  secondary  cyclones  in  group  1 at  the  time  of  cyclo genesis.  The  only 
difference  in  the  application  of  the  grid  to  the  200  mb  chart  was  that  the 
secondary  cyclone  was  centered  in  the  square  marked  X'  (fig.  2). 
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Results 


The  results  of'  Jit  ^tudy  are  presented  in  the  form  of  composite  maps. 

These  maps  dc  r."  t refer  to  a fixed  geographical  region  but  rather  are  con- 
structed with  the  surface  cyclone  center  ns  a common  reference  point.  The 
scale  of  the  map  is  indicated  by  the  2 degree  square  drawn  about  the  refer- 
ence point.  The  area  studied  is  small  enough  so  that  a rectangular  grid  can 
be  used  without  appreciable  distortion  of  the  picture. 

The  composite  200  mb  temperature  distributions  relative  to  the  primary 
surface  cyclone  at  zero-minus-24  hours  are  shown  in  figures  3 and  4 for  the 
group  1 and  group  2 cyclones.  At  this  time  the  differences  between  the  two 
composite  maps  are  rather  slight.  In  both  groups  we  find  the  cold  pool  about 
10  degrees  of  longitude  east  of  the  cyclone  center.  The  warm  pool  appears  to 
the  west  of  the  cyclone  in  group  1 and  west-southwest  of  the  cyclone  in  group  2. 
There  is  as  yet  no  indication  in  figure  3 of  the  protrusion  of  a tongue  of  warm 
air  southeast  of  the  primary  cyclone  or  of  the  formation  of  a secondary  warm 
pool  which  one  might  expect  to  accompany  the  development  of  the  secondary  low 
center  in  th€!  group  1 cyclones. 

At  zero-minus-12  hours  (figs.  5 and  6)  a marked  change  occurs  in  the  tem- 
perature field  of  the  group  2 cyclones.  However,  a glance  at  table  2 shows 
that  this  change  may  be  fictitious  since  only  4 cyclones  were  included  in 
group  2 at  zero-minus-24  hours,  whereas  this  group  includes  seven  cyclones 
at  zero-minus-12  hours. 

A comparison  of  figures  5 and  6 reveals  some  notable  differences  between 
the  group  1 and  group  2 cyclones.  In  the  former  the  cold  pool  appears  about 
19  degrees  of  longitude  to  the  east  of  the  cyclone  while  in  the  latter  the  cold 
pool  is  found  about  the  same  distance  from  the  cyclone,  but  to  the  southeast. 

The  warm  pool  in  group  1 lies  about  5 degrees  of  longitude  to  the  west  of  the 
cyclone.  In  group  2 the  center  of  the  warm  pool  appears  about  10  degrees  of 
longitude  to  the  northwest  of  the  cyclone.  Finally,  there  is  the  obvious  uif- 
ference  in  the  curvature  of  the  isotherms  between  the  two  groups.  The  group  1 
cyclones  show  the  development  of  a tongue  of  warm  air  protruding  soutneast  of 
the  primary  cyclone  from  the  warm  pool.  Tne  group  2 cyclones  show  quite  the 
reverse,  i.e.,  a cold  tongue  extending  northwestward  across  the  cyclone  from 
the  cold  pool. 
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In  retrospect  one  can  now  see  evidence  of  tne  warm  tongue  in  the  secondary- 
producing  lows  as  early  as  zero-minus-24  hours  (fig.  3) • However,  it  does  not 
appear  likely  that  precursory  signs  of  secondary  cyclogenesis  can  be  recognized 
in  the  200  mb  temperature  field  much  more  than  twelve  hours  before  the  dw  eiop- 
ment. 

The  temperature  change  maps  for  the  period  zero-minus-24  hours  to  zero- 
minus-12  hours  (figs.  7 and  8)  also  exhibit  some  interesting  differences  be- 
tween the  group  1 and  2 cyclones.  In  the  former  we  find  warming  at  200  mb 
to  the  east  and  southeast  of  the  primary,  whereas  the  1-tter  show  cooling 
taking  place  in  tne  same  region.  Corresponding  to  the  increasing  temperature 
of  the  warm  pool  in  the  group  2 cyclones,  we  find  rather  considerable  twelve 
hour  warming  nortnwest  of  tire  cyclone  in  figure  3.  This  warming  is  absent  in 
the  group  1 cyclones. 

The  development  in  the  200  temperature  field  associated  with  the  form- 
ation of  a secondary  cyclone  can  be  seen  clearly  in  the  zero  hour  map  (fig.  9)  • 
Here  the  warm  tongue  is  well  developed  southeast  of  the  primary.  The  warm  pool 
is  rotating  about  the  primary  low  and  now  lies  to  the  southwest  of  the  low.  The 
cold  pool  has  continued  to  retreat  eastward  away  from  the  primary  cyclone  c 

The  characteristics  of  the  group  1 cyclones  at  zero  hour  are  even  more 
striking  when  compared  with  those  of  group  2 (fig.  10).  In  group  2 there  is 
no  warm  tongue  development  southeast  of  the  low.  The  warm  pool  has  moved 
southward  relative  to  the  low  but  still  lies  to  the  west-northwest.  The  cold 
pool  shows  no  sign  of  retreating  away  from  the  cyclone. 

The  temperature  change  maps  (figs.  11  and  12)  again  show  marked  differ- 
ences between  the  two  cyclone  groups.  Group  1 cyclones  display  a narrow  band 
of  local  warming  southeast  of  the  primary  cyclone  corresponding  to  the  forma- 
tion of  the  secondary  low  center.  The  warming  in  group  2,  on  the  other  hand, 
is  not  confined  to  one  quadrant  and  appears  to  be  greatest  over  the  surface 
primary  cyclone  itself. 

Twelve  hours  after  the  formation  of  the  secondary  low,  the  temperature 
field  at  200  mb  (fig.  13)  shows  further  development.  The  worm  tongue  east 
and  southeast  of  the  primary  cyclone  becomes  more  pronounced  and  the  warm 
pool  rotates  to  a position  south  of  the  primary.  There  is  some  indication 
of  the  formation  of  a secondary  warm  pool  southeast  of  the  primary  cyclone. 

At  the  same  time  the  group  2 cyclones  (fig.  14)  show  no  further  movement  of 
the  ^arm  pool  relative  to  the  primary  low  and  no  indication  of  the  formation 


Fig.  11.  Average  200  mb  temperature  change  above  group  1 cyclones 
between  zero  - 12  hours  and  zero  hour. 


Fig.  12.  Average  200  mb  temperature  change  above  group  2 cyclones 
between  zero  - 12  hours  and  zero  hour. 


of  a warn)  tongue.  Considerably  more  warming  in  this  period  to  tne  southeast 
of  the  primary  cyclone  is  also  evident  in  the  temperature  change  field  for 
group  1 (fig.  15)  than  for  group  2 (fig.  lb).  However,  figure  15  is  not 
reliable  since  only  three  out  of  the  six  cyclones  were  used  to  construct 
the  zero-plus-12  hour  temperature  field  x'or  group  1. 

The  200  mb  temperature  field  at  zero  hour  relative  to  the  secondary 
cyclones  which  developed  out  of  the  group  1 primary  lows  is  shown  in  figure 
17.  Again  the  warm  tongue  extending  southeastward  from  the  warm  pool  is  evi- 
dent. The  structure  of  the  temperature  field  is  now  not  unlike  that  connect- 
ed with  the  primary  cyclones  in  their  earlier  stages.  The  warm  pool  lies  to 
the  northwest  and  the  cold  pool  to  the  east  or  southeast.  Of  course,  later 
the  warm  pool  will  either  move  to  a position  closer  to  the  new  cyclone  or  a 
new  warm  pool  will  develop  east  of  the  original  one.  This  "jumping"  of  the 
warm  pool  from  the  primary  to  the  secondary  has  been  described  by  Fleagle 
(1947). 

The  composite  maps  were  constructed  from  relatively  few  cases  and  there- 
fore suffer  from  the  usual  deficiencies  of  small  samples.  In  particular  the 
statistical  stability  of  the  averages  is  in  doubt.  Point  values  of  the  stand- 
ard deviation  about  point  means  were  found  to  be  ^s  hign  as  5 degrees  centi- 
grade and  the  average  standard  deviation  was  about  4.5  degrees  centigrade. 

The  large  standard  deviation  indicates  that  no  significance  can  be  attach- 
ed to  small  details  of  the  composite  maps  oi’  to  small  differences  between  the 
group  1 and  group  2 composite  maps.  However,  the  observed  differences  between 
the  group  1 and  group  2 temperature  fields  are  too  large  to  be  ignored,  and  un- 
doubtedly indicate  the  existence  of  different  pnysicai  processes  in  the  two 
groups  of  cyclones. 
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Summary 

The  main  purpose  of  this  study  was  to  determine  the  prognostic  value  of 
the  200  mb  temperature  field  for  tne  prediction  of  secondary  cyclone  forma- 
tion along  the  east  coast  of  the  United  States,  southeast  of  a primary  cyclone. 
The  working  hypothesis  was  that  cyclogenesis  would  be  preceded  by  a charac- 
teristic deformation  of  the  temperature  field  at  200  mb  and  the  appearance 
of  a stratospheric  warm  tongue  or  secondary  war^  pool  southeast  of  the  parent 
cyclone. 

Twenty-four  hours  before  formation  of  the  secondary  it  was  not  possible 
to  distinguish  those  cyclones  which  produced  secondaries  from  those  which  did 
not.  In  both  cases  a 200  mb  warm  pool  was  found  to  the  west  and  a 200  mb  cold 
pool  to  the  east  of  the  primary  surface  cyclcne. 

Beginning  at  twelve  hours  before  formation  of  the  secondary,  the  thermal 
characteristics  of  cyclogenesis  began  to  appear.  A warm  tongue  at  200  mb 
developed  southeast  of  the  primary  low,  the  200  mb  warm  pool  drifted  to  the 
south  of  the  primary  low,  and  the  cold  pool  retreated  toward  the  east.  These 
developments  were  not  observed  in  the  cyclones  which  failed  to  produce  second- 
aries. 

Thus,  although  the  working  hypothesis  appears  to  be  verified,  the  short 
time  lag  between  thermal  developments  at  200  mb  and  surface  cyclogenesis  places 
a practical  limit  on  the  usefulness  of  the  results. 
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VERTICAL  MOTIONS  IN  AN  UPPER  LEVEL  CYCLONE 


Benjamin  M.  Herman  and  John  J.  R.  Kinney 


Introduction 


In  the  last  decade  an  increasingly  dense  network  of  upper-air  observations 
has  become  available  to  the  meteorologist.  These  observations  have  revealed 
the  existence  of  certain  upper-level  phenomena  which  are  not  readily  detect- 
able from  surface  analyses.  One  of  tne  more  important  of  these  phenomena  is 
the  high  level  cold  cyclone.  These  disturbances  are  observed  to  develop  in 
the  southern  end  of  a trough  in  the  westerlies  and,  upon  intensifying,  may 
become  cut-off  from  the  main  stream  and  propagate  down  to  the  surface.  Con- 
sistent feature  of  these  disturbances  are  a pool  of  cold  air  in  the  tropo- 
sphere, a low  tropopause  and  a warm  stratosphere. 

t 

Palman  (1949)  has  described  the  formation  of  a cold  upper-level  cyclone, 
or  cut-off  low,  as  a wave-like  deformation  in  the  upper  flow  which  increases 
in  amplitude  until  the  cold  air  in  the  southern  part  of  the  trough  becomes 
cut  off  from  the  cold  air  to  the  north.  This  cutting  off  may  be  accomplished 
by  subsidence  within  the  cold  air  in  the  central  and  northern  parts  of  the 
trough.  If  the  subsidence  is  intense  enough,  it  may  eventually  result  in  a 
lowering  of  the  top  of  the  cold  air  mass  below  a particular  level,  say  500  mb, 
with  the  cold  air  still  remaining  at  this  level  farther  to  the  south. 

On  upper  air  cheats  the  cut-off  polar  air  appears  as  a cold,  cyclonic 
vortex  in  the  southern  part  of  the  original  upper  trough.  Farther  to  the 
north  a new  belt  of  strong  westerly  winds  can  be  observed.  Thus,  the  whole 
process  can  be  regarded  as  a special  case  of  cyclone  formation  on  the  south 
sid-o  of  the  strongest  westerlies.  A corresponding  transport  of  warm  air 
masses  to  the  north  can  often  be  observed  in  the  formation  of  warm  cut-off 
highs. 

I 

Palmen  and  Nagler  (1949)  have  analyzed  a selected  cold  vortex  of  the 
cut-off  type,  and  have  shown  that,  as  the  cyclone  deepens,  the  temperature 
decreases  within  the  cold  core  in  the  upper  troposphere,  while  there  is  a 
rapid  lowering  of  the  tropopause  (and  subsequently  the  stratosphere)  dir- 
ectly above  the  disturbance.  They  observed  that  a .arm  pool  forms  in  the 
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area  of  Uie  lowering  stratosphere  while  a cold  pool  ionris  in  the  upper  tropo- 
sphere beneath.  It  was  shown  that  the  waim  pocket  could  not  have  teen  a result 
of  any  purely  horizontal  advective  processes,  but  must  have  teen  due  mainly  to 
subsidence.  Similarly,  it  was  concluded  that  the  cooling  of  the  cold  air  in  the 
troposphere  was  brought  about  by  ascending  motion.  (On  the  other  hand,  Nyberg 
(1949)  has  stated  that  in  no  cases  studied  by  him,  did  he  find  any  upward 
motion  of  the  cold  air  in  an  upper  trough.) 

i 

Palroen' s hypotheses  regarding  the  development  of  the  cold  low  and  the 
associated  temperature  field  axe  supported  by  the  analysis  of  Hsieh  (1941)  who 
has  described  not  only  the  development  of  the  low  aloft  but  also  the  downward 
propagation  of  the  upper  level  cyclone  to  the  surface. 

While  numerous  computations  of  vertical  motions  have  been  performed  re- 
cently (see,  e.g.,  Fleagle,  1947),  there  is  as  yet  not  available  a complete 
description  of  the  field  of  vertical  motion  associated  with  the  life  history 
of  an  upper  level  cyclone.  In  this  study  computed  fields  of  vertical  velocity 
at- 500  and  200  mb  are  presented  for  such  a cyclone  during  a two  and  one-half 
day  period  when  the  cyclone  was  moving  southeastward  across  the  north  central 
United  States.  The  500  mb  and  200  mb  levels  were  selected  in  order  to  6how 
the  vertical  motion  fields  in  both  the  troposphere  and  stratosphere. 
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Computai.jon  Methods 

The  method  used  to  compute  vertical  velocities  is  a modified  form  of  the 
graphical  method  of  computing  the  adiabatic  warming  from  constant  pressure 
charts  devised  by  Miller  (1945). 

Miller  has  shown  that 

W * -0.136  T (1) 

where  W is  the  vertical  velocity  of  a particle  with  respect  to  the  isobaric 
surface,  in  cm  sec  , and  ^ T is  the  adiabatic  warming  of  the  air  particle 
in  degree  Centigrade  per  11  hours.  In  order  to  obtain  the  vertical  velocity, 
W,  of  a particle  with  respect  to  the  earth,  the  vertical  velocity,  Wp,  of  the 
isobaric  surface  along  the  horizontal  trajectory  must  be  added  to  W^.  Thus, 

ff  71  TV*  + If  (2) 

t p 

fthen  12  hour  trajectories  are  used,  Wp  can  be  gotten  from  the  equation 

Wp  (cm  sec"1)  = 7 x 10"4  A H(ft)  (3) 

where  AH  is  the  change  in  height  of  the  isobaric  surface  along  the  12  hour 
trajectory.  Therefore  the  vertical  velocity  of  a particle  along  a 12  hour 
trajectory  with  respect  to  the  earth  is  given  by; 

W(cm  sec  1)  = -0.236  AT  + 7xl04  AH 

To  determine  W from  12  hour  traj ectories,  end  points  were  selected  and 
horizontal  trajectories  of  the  sir  particles  were  traced  back  to  determine 
their  origin  12  hours  earlier.  This  was  done  by  moving  the  end  point  back  a 
distance  corresponding  to  6 hours  travel  parallel  to  the  contour  lines  with 
& speed  equal  to  that  of  the  observed  winds  on  the  chart  for  time  t^..  This 
point  was  marked  as  the  mid-point  of  the  12  hour  trajectory  and  was  then 
transferred  to  the  map  for  time  tQ-( t^  - 12  hours).  The  process  was  repeated 
f 07'  another  h hour  period  and  the  final  point  was  marked  as  the  origin  of  the 
12  hour  trajectory  In  areas  where  observed  winds  were  lacking,  gradient 
winds  were  computed. 
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In  and  around  the  low  center,  the  configuration  of  tne  contour  lines  under- 
went considerable  change  during  the  12  hour  periods,  t"us  giving  rise  to  errors 
in  trajectories  determined  in  the  above  n.anner.  In  order  to  minimize  these 

errors,  an  interpolated  mean  chart  was  constructed  for  time  t (t  = t - 6 hours) 

m m i 

and  was  used  as  an  aid  in  determining  trajectories. 

The  t charts  were  constructed  by  averaging  the  temperature  and  height 
m 

values  for  time  t and  t^.  Pilot  balloon  wind  observations  for  time  t were 
o 1 . m 

used  to  supplement  the  interpolated  data.  In  areas  where  no  observed  winds 

were  available  an  average  wind  was  used.  Additional  a ;curacy  was  gained  by 

interpolating  the  t position  of  significant  features  such  as  pressure  centers, 

m 

trough  and  ridge  lines,  and  temperature  centers. 

Twelve-hour  trajectories  in  and  around  the  low  center  could  now  be  com- 
puted with  the  use  of  the  t chart.  The  end  point  was  selected  on  the  t„  chart 

m i 

as  before,  but  moved  back  for  a period  of  only  3 hours  in  the  same  manner  as 

described  above.  T'is  point  was  then  transferred  to  the  t chart  and  the 

m 

particle  moved  back  for  e period  of  6 hours.  This  end  point  was  then  trans- 
ferred to  the  t chart  and  moved  back  the  remaining  3 hours  to  its  origin, 
o 

The  mid-point  of  this  12  hour  trajectory  wet  taken  as  the  mid-point  of  the 
6 hour  trajectory  on  the  t chart. 

lit 

The  value  of  tne  temperature  and  height  at  the  origin  (T  ,H  ) and  end 

o o 

(T^.,Hf)  of  the  12  hour  trajectory  were  noted  and  AH  and  D T were  determined. 

AH,  as  stated  above,  is  the  height  change  of  the  isobaric  surface  along  the 
horizontal  trajectory  and  DT  is  the  isobaric  temperature  change  (T^-T^) 

resulting  from  the  adiabatic  warming. 

The  end  points  of  the  trajectories  were  chosen  at  stations  where  soundings 
were  available,  whenever  possible.  As  this  was  not  possible  for  all  trajectories, 
some  soundings  were  constructed  by  interpolation  from  available  data.  These 
soundings,  therefore,  represent  the  vertical  structure  of  the  atmosphere  through 
the  end  point  of  the  trajectory.  It  was  now  assumed  tnat  the  air-  in  the  layer 
affected  by  the  vertical  motion  moves  with  a constant  horizontal  velocity.  There- 
fore, this  layer  consists  of  the  same  air  particles  at  both  the  beginning  and 
end  points,  thereby  eliminating  advection.  The  temperature  changes  within  this 
layer  are  thus  the  result  of  vertical  motion,  if  ncn-adiabatic  effects  may  be 
considered  negligible. 

To  determine  AT  from  DT,  the  temperature  of  the  beginning  point,  Tq,  was 
plotted  at  the  appropriate  level  (iu  this  study  either  500  mb  or  200  mb)  on  the 
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thermodynamic  diagram,  containing  the  soundings  for  this  trajectory  for  time  t*,» 

X 

The  adiabat  through  T to  the  intersection  with  the  end  point  sounding  was  then 

o 

drawn.  If  the  temperature  at  this  intersection  point  is  denoted  T , then 

AT  =Ta  - Tq  (5) 

A H and  AT  now  having  been  determined,  W,  the  vertical  velocity  of  the 
particle  with  respect  to  the  earth’s  surface  may  be  determined  from  equation 
(4). 

Due  to  the  methods  employed,  the  values  of  IV  thus  obtained  arc  actually 
mean  values  for  the  12  hour  time  interval  and  are  consequently  plotted  at  the 

midpoints  of  the  trajectories.  Also,  since  the  adiabatic  paths  between  T and 

o 

T are  entirely  in  the  layer  immediately  above  or  below  the  level  in  question, 
depending  upon  whether  the  particle  rises  or  sinks,  the  values  of  W represent 
mean  values  for  this  layer,  and  not  actually  for  the  level  in  question.  The 
midpoint  of  this  layer,  however,  is  generally  close  enough  to  the  level  in 
question  so  that  the  deviation  can  be  neglected. 
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The  Vertical  Velocity  Distribution 


On  the  4th  and  5th  of  December  1950  a strong  polar  anticyclone  moved  south- 
eastward from  Canada  into  the  Gret.t-  Plains  of  the  United  States.  At  1230  GMT 
on  the  5th  (fig.  18)  the  polar  anticyclone  was  separated  from  a weaker  ridge  to 
the  east  by  a surface  pressure  trough  in  the  Mississippi  valley,  but  there  was 
no  closed  cyclone  in  the  area.  However,  as  will  be  seen  in  figure  20,  a vex-y 
deep  cyclone  did  exist  aloft  over  the  surface  high.  Forty-eight  hours  later 
(fig.  18)  the  cyclone  was  found  also  at  the  surface  and  the  high  had  split  into 
two  cells,  one  over  the  Gulf  of  Mexico  and  the  cither  over  Quebec.  Subsequently, 
both  the  upper  level  and  surface  cyclones  moved  northward  over  the  Great  Lakes 
while  the  surface  low  deepened  and  the  upper  low  filled  slightly. 

Vertical  velocities  at  tne  500  and  200  mb  levels  in  and  around  tne  upper 
level  cyclone  were  computed  for  the  five  twelve-hour  periods  beginning  at  1500  GMT 
on  the  5th  and  ending  at  0300  GMT  on  the  8th  when  the  cyclone  moved  out  of  the 
U.  S.  radiosonde  network.  The  500  ana  200  mb  maps  and  the  corresponding  vertical 
velocity  distribution  are  shown  in  figures  19  - 29. 

The  solid  lines  on  the  constant  pressure  maps  are  contours  (drawn  for  an 
interval  of  200  geopotential  feet)  and  the  dashed  lines  are  isotherms  (drawn 
for  an  interval  of  4 degrees  Centigrade) . On  the  vertical  velocity  maps  the  thin 
solid  lines  are  isopleths  of  positive  vertical  velocity  and  dashed  lines  are 
ieopleths  of  negative  vertical  velocity  (subsidence).  The  vertical  velocity 
isopleths  (sometimes  called  isanabats)  are  drawn  for  an  interval  of  0.5  cm 
sec~\  Centers  of  rising  motion  are  labelled  R and  centers  of  subsidence  are 
labelled  S.  On  the  200  mb  vertical  velocity  maps  the  tropopause  intersection 
is  shown  as  a heavy  solid  line.  The  interpolated  positions  of  the  cyclones  at 
the  appropriate  isobaric  level  are  indicated  on  the  vertical  velocity  maps  by 
the  letter  L. 

The  vertical  velocity  maps  for  2100  GMT  on  the  5th  (fig.  20)  are  shown  be- 
tween the  constant  pressure  maps  for  1500  GMT  on  tne  5th  (fig.  19)  and  0300  GMT 
on  the  6th  (fig.  21).  During  this  period  the  cyclone,  both  at  200  mb  snd  500  mb, 
moved  southward  and  filled  slightly.  At  500  mb  the  cold  core  of  the  cyclone  be- 
came warmer  while  at  200  mb  the  warm  core  cooled.  However,  the  warm  pool  in  the 
trough  south  of  the  low  center  at  200  mb  warmed  aoout  2 degrees  and  remained 
almost  stationary  over  Kansas. 


Fig.  18.  Surface  weather  maps.  1230  GMT,  5 arid  7 December  1950. 
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The  vertical  velocity  map  for  0900  GMT  on  the  5th  show:.'  three  centers  of 
rising  motion  at  500  mb  (fig.  20)  to  the  east.,  southeast  ar.d  south-southeast 
of  the  lo>v.  West  and  southwest  of  the  low  the  air  was  sinking  at  500  rrb,  while 
in  the  low  itself  the  vertical  velocity  was  almost  zero. 

At  200  mb  (fig,  20)  the  air  was  subsiding  to  the  west  and  southwest  of  the 
cyclone  and  also  to  the  southeast.  However,  in  the  cyclone  itself  the  vertical 
velocity  was  almost  zero. 

Twelve  hours  later,  at  0900  GMT  on  the  6th  (fig.  22),  the  air  at  500  mb 
was  still  rising  to  the  east  and  southeast  of  the  low.  However,  to  the  south 
and  west  of  the  low  subsidence  was  found  at  500  mb. 

The  200  mb  vertical  motion  cnari  for  0900  GMT  on  the  6th  again  shows  sub- 
sidence southwest  oi  the  low  center,  however,  to  the  northwest  the  air  was 
ascending  at  200  mb  over  a region  of  rather  strong  subsidence  at  500  mb.  The 
pattern  of  a stratospheric  lifting  over  tropospheric  sinking  motion  indicates 
the  existence  of  horizontal  convergence  in  the  vicinity  of  the  tropopause, 
which  in  turn  suggests  that  the  low  center  will  move  southeastward  away  from 
the  center  of  high  level  convergence.  A center  of  subsidence  was  also  found 
at  200  mb  over  Tennessee,  well  to  the  southeast  of  the  low.  Below  the  200  mb 
sink  there  existed  ascending  motion  at  500  mb  indicating  the  existence  of 
horizontal  divergence  between  these  two  levels.  This  pattern  again  suggests 
movement  of  the  low  toward  the  southeast,  away  from  the  convergence  and  towards 
the  divergence  aloft.  (A  quantitative  analysis  of  the  divergence  will  be  dis- 
cussed below.) 

Between  0300  and  1500  GMT  on  the  6th  the  low  at  500  mb  (fig.  23)  moved 
southward  with  little  change  in  intensiuy  except  for  an  increase  in  the  gradient 
flow  to  the  east  ana  southeast.  At  200  mb  (fig.  23),  however,  filling  and 
cyclogenesis  occurred  simultaneously.  Although  the  height  of  the  200  mb  sur- 
face generally  increased  over  the  plains,  a new  cyclone  center  appeared  south- 
east of  the  original  low.  The  sinking  motion  of  the  air  at  200  mb  raised  tt-.o 
temperature  of  the  two  pools  of  warm  air.  At  the  same  time  the  gradient  flow 
southeast  of  the  cyclone  increased  considerably  and  the  trough  showed  signs  of 
developing  toward  the  southeast. 

At  2100  GMT  on  the  6th  the  air  was  ascending  at  500  mb  (fig.  24)  in  two 
centers,  one  northeast  and  one  southwest  of  the  cyclone.  The  development  of 
ascending  motion  southwest  of  the  low  is  rather  surprising  since  the  air  had 
been  sinking  in  this  quadrant  at  500  mb.  There  was  also  a marked  change  in 


r'ig.  22.  200  Jib  and  500  mb  vertical  velocity  maps. 

0900  GMT,  6 December  1950. 


I 

Fig.  24.  200  mb  and  500  mb  vertical  velocity  maps. 

2100  GMT,  6 December  1950. 
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the  distribution  of  vertical  motion  at  <i00  mb  at  this  tire  (fig.  24).  Above  tne 
ascending  air  southwest  of  the  low  at  500  mb,  there  appeared  subsidence  in  the 
stratosphere.  A qualitative  consideration  of  the  horizontal  velocity  divergence 
corresponding  to  this  field  of  vertical  motion  indicates  the  existence  of  posi- 
tive divergence  near  the  tropopause  southwest  of  the  low. 

Between  1500  GMT  on  the  6th  and  0300  GMT  on  the  7th  (fig.  25),  the  low  at 
500  "ib  deepened  slightly  and  continued  its  southward  movement.  The  gradient 
steepened  noticeably  in  the  southwest  quadrant  of  the  cyclone.  At  200  mb 
(fig.  25)  the: e was  little  change  in  the  position  or  intensity  of  the  low,  al- 
though the  temperature  in  the  center  increased. 

At  0900  GMT  on  the  7th,  the  center  of  ascending  motion  southwest  of  the 
low  at  500  mb  (fig.  26)  disappeared.  The  dominant  area  of  ascending  motion 
northeast  of  the  low  moved  northward  over  the  Great  Lakes  and  the  upward 
velocities  increased  markedly.  To  the  west  of  the  low  the  air  was  subsiding. 

The  vertical  velocities  at  200  mb  (fig.  26)  were  generally  quite  small  at  this 
time  with  sinking  motion  northwest,  south  and  east  of  the  low  and  rising  motion 
to  the  northeast. 

The  500  mb  low  had  apparently  begun  to  recurve  before  1500  GMT  od  the  7th 
(fig.  27).  At  this  time  the  low  vras  already  moving  northeastward  toward  the 
region  of  ascending  motion  over  the  Great  Lake=>.  Similarly,  the  200  mb  low 
(fig.  27)  began  moving  northeastward.  At  both  levels  there  was  little  change 
in  the  depth  of  the  cyclone. 

By  210*0  GMT  on  the  7th,  the  center  of  subsidence  at  500  mb  (fig.  28),  which 
had  been  west  of  the  low,  moved  to  the  south.  The  center  of  ascending  motion 
moved  eastward  across  the  lakes  and  remained  north-northeast  of  the  cyclone 
center.  At  this  time  (fig.  28)  the  air  at  200  mb  was  sinking  northeast  of  the 
cyclone.  However,  the  vertical  velocities  were  far  greater  at  500  mb  than  at 
200  mb. 

The  low  subsequently  moved  northeastward  both  at  500  mb  and  200  mb,  as 
can  be  seen  on  maps  for  0300  GMT  of  the  8uh  (fig.  29),  and  began  once  again 
to  deepen. 

during  the  period  studied  the  air  at  200  mb  was  generally  sinking  around 
the  low  center.  The  sinking  motion  was  not  confined  to  any  one  quadrant  and 
at  various  times  pockets  of  subsidence  developed  in  all  quadrants  of  the  cyclone. 


Fig.  28.  200  mb  and  500  mb  vertical  velocity  maps. 

2100  GMT,  7 December  1950. 
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At  500  mb  there  appeared  definite  areas  of  rising  as  well  as  sinking  air 
around  the  cyclone.  At  the  beginning  of  the  period,  the  pattern  of  vertical 
notion  at  this  level  resembled  the  well  known  long  wave  tropospheric  model 
with  ascending  motion  to  the  east  and  descending  motion  to  the  wast  of  the 
trough  line.  However  as  the  low  began  to  recurve  northward  toward  tne  end 
of  the  period,  a center  of  rapidly  ascending  air  developed  ahead  (north- 
northeast)  of  the  low  center.  At  the  same  uiiii*;  Hie  uij  to 

rapidly  behind  (south  of)  the  low  center.  By  this  time,  of  course,  the 
north-south  trough  line  had  virtually  disappeared,  the  vortex  was  fully  de- 
veloped, and  the  locus  of  zero  vertical  motion  had  rotated  from  a north- 
south  orientation  to  about  west  northwest  - east  southeast  across  tine  cyclone 
center. 


The  Horizontal  Velocity  Divergence 


According  to  the  equation  of  continuity, 


11  + v_<  f~>  . . 2-  (fw)  = o 

t 


(6) 


where  is  the  air  density,  x,  y and  z are  in  the  east,  north  and  upward 
vertical  directions,  and  u,  v and  w are  the  east,  north  and  upward  vertical 
components  of  the  wind  velocity.  The  continuity  equation  may  also  be  written 
as, 
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If  we  neglect  the  first  terra,  which  represents  the  compressibility  and 
which  is  generally  of  smaller  magnitude  than  the  remaining  terms,  equation 
(7)  reduces  to 

iL-r  = + ^ nv 

^ - — -V  ^ ( 

where  div2W  represents  the  horizontal  divergence. 
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Assuming  a monotonic  change  of  w from  500  mb  to  200  mb,  it  follows  from 
P -,r  (q\  f^at  rising  motion  at  500  mb  and  sinking  motion  at  200  mb  indi- 

cate horizontal  divergence  between  these  levels.  Conversely,  sinking  motion 
at  500  mb  and  rising  motion  at  200  r.b  indicate  horizontal  convergence  between 
the  levels.  In  the  areas  whei e the  vertical  motion  at  both  levels  is  of  the 
same  sign,  the  relative  magnitudes  of  the  vertical  velocities  at  the  two 
levels  determine  the  sign  of  the  divergence. 

Tnere  are  four  possible  combinations  of  vertical  velocities  between  the 
two  levels.  Upon  inspecting  the  vertical  velocity  charts  it  was  noted  that 
the  areas  comprising  these  different  combinations  have  about  the  same  loca- 
tion relative  to  tne  path  of  the  500  mb  low  center  on  each  chart.  In  the 
area  vacated  by  the  low,  200  mb  rising  over  500  mb  sinking  motion  predomin- 
ates. Ahead  of  the  xow  we  find  sinking  at  200  mb  over  rising  at  500  mb. 

To  the  left  of  the  path  (looking  in  the  direction  of  movement)  there  exists 
rising  motion  at  both  levels,  while  to  the  right  of  the  path  we  find  sinking 
motion  at  both  levels.  These  areas  rotate  as  the  path  of  the  low  curves  to 
the  east  and  then  northeast. 

In  order  to  get  some  quantitative  estimate  of  the  magnitude  of  the  di- 
vergence, equation  (8)  was  used  wxth  replaced  by  Au^-,1  . Aw 

was  taken  as  , - w.-_  . ),  and  A z was  taken  as  20,000  ft.  the 

' *00  mb  500  bid  n ' 

average  height  difference  between  the  two  levels.  This  approximation  is 
valid  under  the  assumption  of  a lineal1  change  of  w with  height.  Computa- 
tions were  made  at  points  of  maximum  difference  of  vertical  velocity  between 
the  two  levels.  These  points  along  with  the  magnitudes  of  div,  \s/  in  sec  ^ 
are  shown  on  the  divergence  charts  (figs,  30  - 34).  Thin'  divergence  may  be 
interpreted  as  the  mean  divergence  in  the  layer  between  500  and  200  mb. 

In  general  it  can  be  seen  from  the  divergence  charts  that  the  low  moves 
away  from  an  area  of  convergence  and  toward  an  area  of  divergence.  (The  six- 
hourly  cyclone  positions  shown  on  the  divergence  maps  are  for  the  o00  mb  level) . 
An  interesting  example  of  this  can  be  seen  on  the  map  for  0900  GUT  of  the  7th 
(fig.  33).  The  area  of  convergence,  previously  to  the  nortn  of  tne  low  center 
(figs.  30,  31,  32)  moved  to  the  west  and  southwest  of  tne  low  along  the  path 
previously  taken  by  the  low.  This  was  reflected  in  a change  in  the  path  of 
the  low  towards  the  southeast  and  east.  On  the  folic. ung  chart  (fig.  34)  It 
can  be  seen  that  the  area  of  convergence  intensified  and  moved  to  the  southeast. 
The  low  then  recurved  towards  the  northeast  and  finally  to  the  north,  directly 
into  an  area  of  strong  divergence  over  the 


Great  Lakes 


0900  GicT,  t>  December  1950. 


Q O' 


Fig.  32.  Horizontal  divergence  between  200  and  500  mb. 
2100  GMT,  6 December  1950. 


Fig.  33.  Horizontal  ‘divergence  between  200  and  500  mb. 
0900  GOT,  7 December  1950. 
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COMPUTATIONS  OF  LONG  WAVS  KADIAflONAL 
COOLING  IN  EXTRATROPICAL  CYCLONES 


GEORGE  OUR.  NG 


Introduction 

There  has  been  much  research  on  mean  infrared  radiational  cooling  In  the 
atmosphere.  These  studies,  wnich  present  an  average  picture  of  atmospheric 
radiational  processes,  tend  to  smooth  out  the  irregularities  which  occur 
daily  in  various  synoptic  situations.  One  of  the  most  important  synoptic  sit- 
uations in  need  of  further  understanding  is  the  cyclone.  The  aim  of  'inis  paper 
is  to  determine  tne  effect  of  the  synoptic  conditions  within  a storm  upon  the 
radiational  cooling  rates.  This  in  turn  may  lead  to  a better  understanding 
of  the  role  of  radiational  cooling  on  cyclone  structure  and  development. 


Summary  of  Earlier  Work 

The  greater  part  of  previous  research  on  infrared  radiational  cooling  has 
been  concentrated  upon  the  computation  of  cooling  rates  for  mean  or  idealized 
conditions.  For  clear  Arctic  conditions,  the  mean  summer  maximum  cooling  rate 
is  1.2°c/day  and  the  mean  winter  maximum  cooling  rate  is  0.8°C/day  (Penner,1948) . 
For  middle  latitudes,  average  clear  sky  cooling  rates  have  a maximum  of  1°  to 
2 C/day  (Moller,  1950  and  Lcuion,  1952).  In  clear  equatorial  air,  mean  maximum 
cooling  rates  are  2 to  3°C/day  (Elsasser,  1940). 

Studies  of  radiational  cooling,  assuming  mean  cloud  conditions  as  well  as 
mean  lapse  rates,  indicate  that  with  the  presence  of  clouds  the  cooling  rate 
increases.  In  an  average  middle  latitude  cyclone  the  cooling  rate  at  tue  top 
of  tne  lower  cloud  layer  is  2°C/day  and  at  tne  top  of  the  higher  overcast  is 
5°C/day  (Moller,  1950) . V'ritk  mean  cloudy  conditions  in  the  Nortuern  Hemisphere 
the  cooling  rate  is  greater  than  2° c/day  at  the  tops  of  the  clouds  (London, 1952) . 

The  theoretical  infrared  radiational  properties  of  a cloud  layer  50m  thick. 

It 

in  a normal  atmosphere,  were  investigated  by  Roller  (1943)  • At  a height  of  2 km 


»•  w kA 
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the  cloud  top  cooled  at  a rate  of  to. 5°  c/hr  and  the  base  warmed  at.  a rate  of 
1.7°C/hr.  fthen  placed  at  6 km.  the  cloud  top  cooled  at  a rate  of  9.6°C/hr 
and  the  base  warmed  at  a rate  of  to. 9 C/hr.  Moller  noted  that  suen  a large 
cooling  at  the  top  and  a warming  at  the  base  would  result  in  a very  unstable 
lapse  rate  within  the  cloud. 

Synoptic  studies  of  long  wave  radiational  cooling  have  been  made  by  Alder 
(1951)  and  by  Neiburger  and  Thompson  (1953).  Alder  investigated  a polar  out- 
break over  the  United  States  and  reported  that  maximum  cooling  rates  were  found 
above  clouds,  that  maximum  heating  in  the  vicinity  of  clouds  occurred  at  the 
bases,  that  atmospheric  heeting  was  found  at  the  bases  of  inversions,  and  that 
atmospheric  cooling  was  found  at  the  tops  of  inversions.  In  their  study  of 
the  subtropical  anticyclone  in  bhe  Pacific  Ocean,  Neiburger  and  Thompson  found 
maximum  cooling  rates  greater  than  2-5°C/day  at  the  700  mb  level.  These 
maximum  rates  cccu.  red  at  the  top  oi  a subsidence  inversion  in  the  eastern 
part  of  the  high  pressure  cell. 

It  is  known  that  radiational  cooling  in  the  free  air  depends  upon  the  tern- 
peratoire  distribution,  the  moisture  distribution  and  the  cloudiness  distribu- 
tion. In  a cyclonic  situation  the  most  important  of  these  three  is  in  the  ef- 
fect of  cloudiness.  From  previous  studies  of  longwave  radiational  cooling  and 
theoretical  reasoning  one  would  expect  to  find  large  cooling  at  the  top  surfaces 
of  the  cloud  layers  within  the  cyclone.  Also,  one  would  expect  to  find  a general 
warming  at  the  bases  of  the  cloud  layers.  At  the  frontal  inversions  in  a cyclone 
the  temperature  distribution  influences  the  cooling  rates  to  a great  extent.  The 
rate  of  temperature  change  would  be  expected  to  change  from  a warming  rate  at  the 
base  of  the  inversion  to  a cooling  rate  at  the  top  of  the  inversion. 
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Procedures  and  Methods  of  Computation 


In  all  methods  of  computing  radiational  temperature  changes,  the  cooling 
is  determined  from  the  divergence  of  the  net  radiational  flux.  The  equation 
for  the  cooling  is  given  by 

AX  _ 

Af  (1) 


where  AT/  At  is  the  local  temperature  change  per  unit  time 
g is  the  acceleration  of  gravity 


c is  the  specific  heat  of  dry  air  at  constant 
P 

pressure 

ASn/Ap  is  the  change  of  net  radiational  flux  S witn 
cnanging  pressure  in  the  vertical  direction. 


This  formula  is  derived  from  the  first  law  of  thermodynamics  assuming  that  the 
effects  of  compression  and  advection  can  be  neglected. 

There  are  direct  and  indirect  metnods  for  computing  radiational  cooling. 
Eisasser's  method  (1942),  for  exam; le,  is  an  indirect  one  and  requires  a deter- 
mination of  the  net  flux  graphically  and  a further  computation  to  obtain  the 
cooling.  This  method  is  not  very  practical  for  synoptic  applications  since 
it  tends  to  smootn  out  irregularities.  Bruinenberg' 6 tabular  method  (1946), 
which  is  based  upon  a simplified  water  vapor  spectrum,  obtains  the  cooling 
directly  from  an  integration  of  tne  flux  divergence  above  and  below  the  level 
at  which  the  cooling  is  being  computed.  In  this  study,  D.  L.  Brooks'  tabular 
method  (1950)  was  used,  ratner  than  Bruinenberg' s,  because  of  its  availability. 
Brooks  derived  his  method  from  Bruinenberg1 s equation  (1946)  lor  the  instan- 
taneous rate  of  temperature  change  due  to  infrared  radiation.  The  essential 
difference  between  the  two  methods  is  that  Brooks  makes  . ;e  of  empirical  emis- 
sivity  measurements  rather  than  the  simplified  water  vapor  spectrum  theorized 
by  Bruinenberg. 
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Brooks’  working  equation  for  the  instantaneous  rate  of  temperature  cnange 
due  to  water  vapor  radiation  at  a reference  level  in  the  atmosphere  is 


37- 

<?fc 


A 


1 - i-f**ry  (.p j r 

L Pur  ' ^4;  - v 


where  is  the  local  temperature  change  in  degrees  centigrade  pe; 

minute  at  the  reference  level 


q is  the  specific  humidity  at  the  reference  level 


c is  the  specific  heat  of  dry  air  at  constant  pressure 


p 2/2 

( /1000)  ' is  the  square  root  pressure  correction  factor  to 

obtain  the  corrected  optical  path  length  at  the 
reference  level  where  the  pressure  is  p.  The  need 
for  a pressure  correction  factor  is  explained  in 
the  section  on  discussion  of  assumptions. 


( 6 T ) is  the  black  body  radiation  of  a slab.  ( is  the 

Stefan  - Boltzmann  constant.) 


S is  the  infrared  radiational  flux  . 


£ is  the  emissivity  of  water  vapor  &uu  is  defined 
as  the  actual  emission  divided  by  the  black  body 
emission.  It  is  a function  of  the  corrected 
optical  path  length  w . 

» is  the  corrected  optical  path  length  for  water 
vapor  . 

w is  the  value  of  cj  s/c at  the  upper  boundary  of 
the  radiating  atmosphere. 


In  equation  (2),  the  summation  over  n is  for  layers  above  the  reference 
level  into  which  the  atmosphere  has  been  divided,  and  the  summation  over  m 
is  for  layers  below  the  reference  level.  The  points  on  the  radiosonde  sound- 
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ing  are  considered  to  bound  the  layers  into  which  the  atmosphere  nas  been 
divided  for  puiposes  of  the  summation  of  the  individual  contributions.  Mtnin 
tiie  brackets  of  equation  ( <) , A ( <5>  T^)  is  tiie  change  in  ( 6 T^)  in 

each  1,-yer,  3 ' /of*  is  tiie  average  slope  of  the  emissivity  curve  for  an  incre- 
ment  of  w for  eacn  layer,  and  o)i^/o)w  is  equal  to  ( 6"  T**)  (c)  £.  /c)w  evalu- 

ated at  the  top  of  the  water  vapor  atmosphere.  The  term  cjS,-./o)w  becomes 
zero  at  a black  body  surface  such  as  the  ground  or  a cloud  sine?  w becomes 
infinite  and  d)<£  f c}w  becomes  zero. 

This  method  was  applied  in  the  following  manner.  The  first  step,  as  in 
all  methods  of  computing  radiative  fluxes  and  cooling  rates,  is  to  determine 
w,  the  corrected  optical  path  length,  for  each  layer  of  the  sounding.  For 
small  layers  this  is  given  by  (Brooks,  1950) 


* rented)  = - q ^ p ( A 1/2  (3) 

g ^ 1000  / 


where  w is  the  cerrpnt.fid  optical  path  length  in  g/ cm 


q is  the  mean  specific  humidity  for  tiie  layer  in  g/kg 
p is  the  thickness  of  the  layer  in  mb 

2 

g is  the  acceleration  of  gravity  in  cm/sec 
p is  the  mean  pressure  of  the  layer. 


Next,  the  increments  of  path  length  (w)  are  summed  up  with  respect  to 
the  level  at  which  the  temperature  change  is  being  computed,  i.e„,  the  refer- 
ence level.  At  the  reference  level  w is  taken  as  zero  and  increases  in  either 
direction.  From  Brooks'  tables  tiie  values  of  A ( <£>'  T^) , (d)  £ 75 *>  and 
o)J|ycJw  are  obtained.  If  there  are  no  cloud  decks,  the  summation  of  equation 
(2)  extends  to  the  top  of  the  water  vapor  atmosphere.  It  is  sufficient  to  go 
up  as  high  as  200  mb  or  150  mb  since  above  these  levels  there  is  very  little 
contribution  to  the  cooling  at  lower  heights.  If  there  is  a cloud  deck  the 
summation  is  carried  out  up  to  the  base  of  the  cloud.  Since  the  cloud  is  a 
black  body,  the  term  c)£^/o)w  becomes  zero  and  does  not  have  to  be  included 
in  equation  (2). 


The  summation  of  equation  (2)  is  then  carried  out  and  the  rate  of  tem- 
perature change  is  obtained.  For  partially  cloudy  skies  a weighted  cooling 


rate  was  obtained  by  weighting  the  cooling  rates  computed  for  clear  skies  and 
total  cloudiness  according  to  the  amount  of  clouds  present.  Tnis  weighted 
cooling  rate  formula  is 
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cloudy 


clear 


(4) 


where  <)  T/o)  t is  the  weighted  cooling  rate, 


N is  the  amount  of  clouds  in  tenths, 


<■*  T/<*Hloudy 


is  the  cooling  rate  assuming  total  cloudiness. 


and 


(o)  TAH)clear  the  cooling  rate  assuming  clear  skies. 


L "5  r-xdiitior.nl  cooli..6  ioi.es  minin  tne  clouds  were  computed  in  the 

following  manner.  Assuming  thft  each  50  mb  layer  within  a cloud  acts  as  a per- 
fect black  body,  the  cooling  Is  given  by 

AT  = K A 

At  A P (5) 

is  the  rate  of  temperature  change  per  day,  and 

is  the  curvature  of  the  temperature  sounding  . 

is  a constant  which  depends  upon  the  thickness  of  cloud  chosen 
as  a complete  black  body  radiator  (in  these  calculations,  50  mb), 
and  upon  the  mean  temperature  of  the  layer. 

The  method  by  which  equation  (5)  can  be  derived  is  discussed  by  Brunt  (1952). 

The  thickness  of  cloud  that  will  behave  as  a complete  black  body  is  not  pre- 
cisely known  at  the  present  time.  However,  no  matter  what  thickness  is  assum- 
ed, it  will  only  change  the  magnitude  of  the  radiational  cooling,  not  the 
sign  or  the  distribution. 


A T 

where  t 


-1—2 
A P 
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From  the  previous  discussion  of  the  method  for  obtaining  the  radiational 
cooling,  it  can  be  seen  that  three  quantities  must  be  known  in  order  to  carry 
out  the  computation.  These  are 

1)  The  temperature  distribution  in  the  vertical 

2)  The  moisture  distribution  in  the  vertical 

3)  The  amount  of  clouds  and  heights  of  bases  and 
tops  of  clouds. 

Temperatures  were  obtained  from  raiio sonde  data  and  plotted  on  tephigrair.s. 
Where  the  temperature  soundings  did  not  extend  up  to  li>0  mb  they  were  extra- 
polated to  tnat  point  on  the  basis  of  the  temperature  fields  on  the  200  mb  and 
150  mb  facsimile  maps.  Dew  points,  obtained  from  radiosonde  data,  were  used 
as  a measure  of  the  moisture  distribution.  They  were  extrapolated  to  150  mb 
assuming  a humid  stratosphere  with,  saturation  at.  the  tropopause  (see  Hoffer, 
1951) . All  the  assumptions  made  in  this  section  are  discussed  below. 

Cloud  distributions  were  based  on  synoptic  data,  three-hourly  data,  pilot 
reports  and  radiosonde  data.  Wherever  cloud  data  were  lacking,  cloud  distri- 
butions were  determined  by  using  the  soundings  and  on  the  basis  of  continuity 
in  a cyclonic  cloud  distribution. 
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Discussion  of  Assumptions 

The  warming  of  the  atmosphere  due  to  absorption  of  short  wave  solar  radia- 
tion was  ignored  in  this  study  because,  for  the  most  part,  it  is  negligible 
compared  to  the  temperature  changes  due  to  long  wave  radiation.  In  a clear 
atmosphere,  wanning  by  absorption  of  solar  radiation  amounts  to  only  a few 
tenths  of  a degree  per  day.  With  cloudy  skies,  according  to  a sample  calcula- 
tion, the  absorption  of  solar  radiation  may  produce  a warming  of  2.5°  C/day 
in  the  top  20  mb  of  the  cloud.  This  warming  turns  out  to  be  almost  as  large 
as  the  infrared  cooling  rates  computed  in  this  paper  near  the  cloud  tops. 
However,  it  should  be  pointed  out  that  the  long  wave  cooling  rates  at  the 
top  of  the  cloud  are  much  larger  than  those  presented  in  this  paper,  just 
as  the  temperature  of  the  ground  surface  is  muen  higher  tnan  the  temperature 
of  the  air  a few  centimeters  above  it  during  the  daytime.  Also,  solar  radia- 
tional  processes  operate  only  while  the  sun  is  shining,  while  long  wave  radia- 
tional  processes  operate  througnout  the  day. 

The  extrapolation  of  temperature  on  the  basis  of  the  temperature  field  on 
upper  level  constant  pressure  charts  is  fairly  accurate  since  the  temperature 
field  is  continuous.  The  extrapolation  of  dew  points  up  to  150  mb  is  based  upon 
the  presumption  of  a humid  stratosphere.  At  the  present  time  our  knowledge  of 
the  stratosphere  leads  us  to  believe  that  there  is  an  equal  probability  of  there 
being  a humid  or  dry  stratosphere.  Hoffer  (1951)  computed  tropospheric  radia- 
tional  cooling  curves  for  both  dry  and  humid  stratospheric  conditions.  His  re- 
sults indicate  that  there  is  little  difference  between  the  tropospheric  cool- 
ing computed  from  both  assumptions.  The  humid  stratospheric  assumption  was 
chosen  in  this  study  as  a matter  of  mechanicaJL  convenience  in  the  extrapolation 
of  dew  points. 

a pressure  correction  factor  is  necessary  in  the  atmosphere  because  it  has 
been  found  that  the  absorption  coefficients  of  water  vapor  deperd,  in  a rather 
complex  manner,  upon  the  air  pressure.  Previous  investigators  have  used  a square 
root  pressure  correction  factor  (sec,  for  example,  E.lsasser,  19A2) . This  type 
of  pressure  correction  was  used  in  this  paper.  Recent  studies  (see,  for  example, 
Kaplan,  1952  ) indicate  that  a linear  pressure  correction,  more  correctly 
defines  the  pressure  uependence  of  absorption.  This  new  correction  would  not 
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influence  the  results  of  this  study  except  at  high  levels  and  even  there  it 
would  not  change  the  pattern  of  radiative  cooling. 

The  long  wave  radiational  cooling  of  carbon  dioxide  in  the  atmosphere 
was  ignored  in  this  study  because  it  is  insignificant  when  compared  to  the 
water  vapor  radiational  cooling.  According  to  Brooke  (1950),  CO  2 cooling 
is  generally'  about  2>%  as  large  as  water  vapor  cooling. 

Some  of  the  compute  tiros  made  in  this  study  with  Brooks’  method  were 
checked  by  the  use  of  Bruinenberg’ s tabular  method.  The  largest  difference 
in  cooling  rates  computed  by  both  methods  was  about  0.7°C/day.  When  Brooks 
(1950)  compared  his  method  with  Bruinenberg* s,  he  found  that  bruin en berg’ s 
tables  gave  consistently'  lower  values  — generally  about  2/3  as  great.  The 
reason  for  this  difference,  according  to  Brooks,  is  the  variation  in  emissiv- 
ity  with  temperature  derived  theoretically  by  both  of  these  researchers.  Our 
present  knowledge  cannot  tell  us  which  of  these  methods  gives  more  accurate 
results. 


gyp optic  Situations  Investigated  and  Sources  of  Data 

Cooling  rates  were  computed  in  two  different  cyclones.  The  first  cyclone 
(No.  I)  was  centered  in  southern  Illinois  at  0330  GMT,  1 April  1953.  Within  a 
period  of  24  hours  it  moved  eastward  across  the  country  to  northeastern  Virginia 
and  occluded  (figs.  35,  36,  37).  In  general,  there  were  overcast  sky  conditions 
and  rain  areas  north  of  the  cyclone  center.  Cooling  rates  were  computed  above 
Buffalo,  New  iork;  Mount  Clemens,  Michigan;  Rantoul,  Illinois;  and  Omaha,  Nebraska, 
for  0300  GMT,  1 April  1953.  For  1500  GMT,  cooling  rates  were  computed  above 
Rantoul,  Mount  Clemens,  Buffalo,  Pittsburgh,  Pennsylvania  and  Washington,  D.C. 

For  0300  GMT  2 April  1953,  cooling  rates  were  computed  above  Dayton,  Ohio, 
Pittsburgh,  and  Lakehurst,  New  Jersey. 

The  second  storm  (No.  II)  was  centered  in  eastern  Kansas  at  0300  GMT, 

20  November  1953.  In  12  hours  it  moved  northward  to  northeastern  Nebraska  and 
its  central  pressure  dropped  from  1002  mb  to  996  mb  (figs.  38,  39).  Widespread 
cloudiness,  rain  and  snow  were  prevalent  at  the  center  f the  storm  and  to  the 
north  of  the  center.  For  0300  GMT,  20  November  1953,  cooling  rates  were  computed 


above  Oklahoma  City,  Oklahoma;  Columbia,  Missouri;  and  Rantoul.  for  1500  GMT 
cooling  rates  were  computed  above  International  Falls,  Minnesota;  Saint  Cloud, 
Minnesota;  Omaha  ar.d  Columbia. 


Results 

The  most  important  effect  on  the  radiational  cooling  in  a cyclone  is  caused 
by  the  dir t”i bution  of  cloudiness  within  the  storm.  With  cloudy  skies,  the  air 
below  tne  clouds  does  not  cool  as  rapidly  as  it  would  with  clear  skies.  The 
air  immedittely  above  the  clouds  cools  at  a greater  rate  than  it  would  with 
clear  skies.  These  two  effects,  slight  cooling  or  warming  immediately  below 
the  clouds  in  a cyclone,  and  large  cooling  immediately  above  the  clouds,  can 
be  seen  in  the  cooling  rate  curves  (figs.  40  - 44). 

The  cooling  rate  curves  shew  the  temperature  and  dewpoint  soundings  ana 
the  computed  cooling  rates  on  one  diagram.  The  positions  and  amounts  of  clouds, 
and  the  frontal  zones  are  also  indicated  on  these  diagrams.  Wherever  there  is 
partial  cloudiness  the  cooling  rates  are  weighted  cooling  rates. 

Large  cooling  rates  are  found  at  the  tops  of  the  overcast  layers.  In 
cyclone  I (figs.  40  - 42)  cooling  rates  of  3 C/day  are  found  above  the  thick 
altostratus  - altocumulus  overcasts  of  the  storm.  Warming  rates  in  cyclone  I 
(figs.  40  - 42)  are  equal  to  about  1°  C/day  at  the  bases  of  the  cloud  decks  of 
the  storm. 

The  cooling  curves  for  cyclone  II  are  shown  in  figures  43  and  44.  The 
most  striking  feature  o?  these  curves  is  the  cooling  rate  of  more  than  8°  c/day 
at  the  top  of  the  cloud  layer  at  Columbia  (fig.  43).  This  large  rate  resulted 
from  the  fact  that  there  was  a steep  inversion  at  the  top  of  the  cloud.  The 
cooling  rate  was  computed  at  the  top  of  this  inversion  which  was  about  15  mb 
above  the  cloud  top.  This  level  is  at  a higher  temperature  than  the  cloud 
top  or  any  ^oiut  above  the  inversion  and  will  therefore  lose  a great  deal  of 
heat  by  radiation.  Cooling  rates  at  the  tops  of  the  cloud  systems  in  this 
cyclone  are  also  generally  about  3°  C/day.  Warming  of  as  much  as  1.8°  C/day 
(fig.  44)  is  found  immediately  below  the  cloud  bases  at  Omaha. 
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Fig.  40.  Lapse  rate  and  cooling  rate  curves.  0300  GMT,  1 April  1953. 
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Fig.  41.  Lapse  rate  and  cooling  rate  curves.  1500  GMT,  1 April  1953. 


FiK>  43*  Lapse  rate  and  cooling  rate  curves.  0300  GMT,  20  November  1953. 


Lapse  rate  and  cooling  rate  curves.  1500  GMT,  20  November  1953 
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The  effect  of  frontal  inversions,  which  cause  rediational  cooling  at  their 
tops  and  warming  at  their  bases,  can  also  be  seen  in  the  cooling  areas.  Ex- 
amples of  these  large  changes  from  warming  to  cooling  £3  car  from  the  base 

of  the  frontal  inversion  to  the  top  can  be  seen  at  Mount  Clemens  in  figure  4C 
and  at  Oklahoma  City  in  figure  43. 

Within  clouds  cf  normal  temperature  distribution,  there  is  warming  near 
the  base,  cooling  near  the  top  and  a tendency  for  zero  cooling  throughout  the 
remainder  of  the  cloud. 

Irregular  temperature  soundings  caused  the  erratic  cooling  and  wanning 
rates  above  Buffalo  (fig.  40). 

As  was  indicated  previously  the  radiational  cooling  is  a function  of  the 
temperature,  the  moisture,  and  the  cloud  distributions.  The  effects  of  cyclonic 
cloudiness  on  the  cooling  rates  wore  discussed  in  the  previous  section.  Now 
the  effects  of  temperature  and  moisture  variations  upon  the  radiational  cooling 
rates  shall  be  discussed. 

It  can  be  6hown  that  the  radiational  cooling  depends  upon  the  change  of  the 
temperature  lapse  rate.  It  is  to  be  expected,  therefore,  that  maxima  and  minima 
of  cooling  would  occur  at  places  of  large  curvature  of  the  temperature  profile, 
such  as  inversions.  Large  cooling  would  be  expected  at  the  tops  of  inversions 
and  wanning  would  be  expected  at  the  bases  of  inversions.  These  effects  have 
already  been  shown  fox-  frontal  inversions.  This  temperature  dependence  of 
radiational  cooling  can  be  seen  in  the  cooling  curves. 

Figure  40  contains  the  Buffalo  cooling  curve  for  0300  GMT,  1 April  1953. 
Maximum  coding  rates  are  found  at  the  tops  of  all  four  inversions  in  the  tropo- 
sphere and  maximum  warming  rates  are  found  at  the  bases  of  these  inversions. 

Figure  43  contains  the  cooling  curves  for  itantoul  at  0300  GMT,  20  November 
1953.  Through  the  steep  inversion  from  850  mb  to  840  mb  the  rate  of  temperature 
change  varies  from  a warming  rate  of  4°  C/day  at  the  base  to  a cooling  rate  of 
over  6”  G/day  at  the  top.  A 5°  C/day  cooling  rate  at  948  mb  is  due  to  the  fact 
that  this  point  is  at  the  top  of  a ground  inversion. 

The  dependence  of  radiational  cooling  upon  moisture  content  can  also  be  seen 
in  the  cooling  curves.  A moist  layer  radiating  upward  to  a dry  atmosphere  will 
cool  substantially.  A dry  layer  radiating  to  a moist  atmosphere  will  not  cool  as 
effectively.  In  figure  42,  in  the  layers  above  the  cloud,  the  curvature  cf  the 
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temperature  profile  at  Dayton  does  not  change  substantially  and  therefore  does 
not  change  the  amount  of  radiational  coding.  At  615  mb,  the  top  of  a moist 
layer,  the  cooling  rate  is  3-7°  c/day.  At  595  mb,  the  top  of  the  dry  layer 
that  is  located  immediately  above  the  moist  layer,  the  cooling  rate  decreases 
to  1.9°  C/day.  This  change  in  cooling  rate  is  due  to  the  sharp  change  in  the 
dew  point  lapse  rate  in  the  layers  involved. 


Conclusions 

Just  above  overcast  cloud  layers  the  cooling  is  quite  large,  generally 
2 to  4°  C/day.  Cloud  tops  behave  as  black  bodies  and  radiate  very  efficient- 
ly at  their  own  temperatures.  The  downwexd  radiation  from  the  water  vapor 
above  the  cloud  is  not  nearly  enough  to  compensate  for  the  upward  radiation 
of  the  cloud  top.  This  results  in  a larg6  net  flux  directed  upward  at  the 
cloud  top.  This  net  flux  increases  rapidly  with  height  immediately  above  the 
cloud  top,  since  the  upward  flux  decreases  at  a lower  rate  than  the  downward 
flux.  These  processes  cause  a large  positive  flux  divergence  immediately 
above  the  Gloud  and  therefore  a large  cooling  rate. 

Immediately  below  cloud  decks  there  is  slight  radiational  heating,  about 
1°  C/day.  There  is  heating  in  the  vicinity  of  cloud  bases  because  the  diver- 
gence of  the  net  flux  is  negative,  i.e.,  there  is  a convergence  of  the  net 
flux  immediately  below  the  cloud  base.  This  convergence  of  net  flux  results 
from  the  fret  that  the  upward  radiation  from  the  layers  immediately  below  the 
cloud  decreases  more  rapidly  with  increasing  height,  than  does  the  downward 
radiation  from  the  cloud  and  layers  immediately  underneath  it.  It  can  easily 
be  seen  that  radiational  processes  tend  to  make  the  air  within  the  clouds  more 
unstable,  perhaps  aiding  in  the  further  development  of  the  cloud. 

There  is  a tendency  for  aero  radiational  cooling  to  take  place  in  the 
interior  of  clouds.  This  neutral  condition  is  caused  ty  the  fairly  constant 
lapse  rates  generally  observed  within  clouds.  With  constant  lapse  rates,  there 
is  no  curvature  of  the  temperature  profile  and  therefore  no  cooling. 
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Fronts,  inversions  and  discontinuities  of  lapse  rate  in  general  are  charac- 
terized by  a change  in  the  curvature  of  the  temperature  profile.  This  results 
in  great  changes  in  the  cooling,  since  the  radiational  cooling  is  a function 
of  the  curvature  of  the  temperature  profile.  Large  cooling  rates  are  found 
at  the  tops  of  frontal  inversions  and  other  inversions.  Warming  is  generally 
found  at  the  bases  of  these  inversions.  The  magnitudes  of  these  rates  of  tem- 
perature change  depend  upon  the  steepne, \ of  the  inversion.  In  this  study,  the 
extremes  ranged  from  +4°  C/ day  <it  the  base  of  an  inversion  to  -6°  c/day  at  the 
top  of  an  inversion.  Thus  long  wave  radiational.  processes  tend  to  destabilize 
the  air  or  destroy  the  inversion.  However,  the  magnitude  of  the  radiational 
destabilizing  process  will  decrease  in  intensity  as  the  inversion  is  destabilized. 

The  large  radiational  cooling  rates  found  in  this  study,  especially  at  cloud 
tops  and  at  the  tops  of  inversion  layers,  may  help  to  explain  some  of  the  erratic 
temperature  changes  which  are  occasionally  observed  on  upper  air  maps.  They  also 
suggest  a possible  source  of  error  in  the  computation  of  vertical  velocities  by 
the  adiabatic  method  as  well  as  in  the  advective  temperature  forecasts. 

It  is  clear  that  the  role  of  radiational  cooling  is  to  "smooth  out"  the 
lapse  rate  curve.  Inversions  tend  to  be  destroyed  by  the  large  radiational  cool- 
ing at  the  top  and  the  radiational  warming  at  the  bottom.  Thus,  soundings  on 
which  appear  a multiplicity  of  inversion  layers  are  probably  quite  transitory 
since  they  are  const&Oly  undergoing  modification  by  radiational  processes. 

The  alteration  of  inversion  layers  by  radiational  cooling,  however,  is 
greatly  modified  by  the  existence  of  clouds.  A cloud  layer  with  its  top  at  the 
base  of  the  inversion  layer  (e.g.,  stratus  and  stra to cumulus)  cause  radiational 
cooling  to  prevail  at  the  base  of  the  inversion  and  so  maintain  the  inversion 
layer.  Similarly,  a cloud  layer  with  its  base  at  the  top  of  an  inversion  layer 
(e.g.,  frontal  altostratus)  may  prevent  the  destruction  of  the  inversion  by 
diminishing  the  radiative  cooling  which  generally  prevails  at  tops  of  inversions. 

Finally,  the  cooling  curves  indicate  that  inversions,  including  frontal  in- 
versions, may  move  up  or  down  in  response  to  radiative  cooling.  This  fact  may 
explain  the  difficulty  which  analysts  often  experience  in  establishing  continuity 
between  vertical  crosB  sections.  Apparently  frontal  surfaces  in  the  free  atmoshere 
are  not  always  substantial  surfaces. 
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MONTHLY  FREQUENCIES  OF  CZCLOGENESIS  IN  THE  EAST 
COASTAL  REGION  OF  THE  UNITED  STATES 


Monthly  frequencies  of  east  coastal  cyclogenesis  have  beeu  determined  by 
Miller  (1946)  for  the  months  October  through  April  on  the  basis  of  data  from 
the  period  1929-1939*  In  the  present  study  MiJLLer' s data  have  been  supple- 
mented by  the  addition  of  cyclones  from  the  ten-year  period  1942-1952.  The 
cyclone  season  has  also  been  extended  to  include  the  months  September  and  May. 

For  the  purpose  of  this  study,  the  east  coastal  region  of  the  United 
States  was  defined  by  the  same  boundaries  employed  by  Miller  (see  Miller's 
figure  1) . 

The  cyclone  types  described  by  Miller  were  adopted  in  this  study.  These 
include  the  type  A or  cold  frontal  cyclone  and  the  type  B or  warm  frontal 
cyclone,  both  of  which  are  illustrated  in  Miller's  paper.  It  was  found  necess- 
ary, however,  to  Introduce  another  cyclone  type,  called  type  C , which  forms  at 
the  peak  of  the  warm  sector  (or  point  of  occlusion)  of  an  occluded  cyclone. 

Type  C cyclones  may  also  form  on  or  near  north-south  oriented  cold  fronts. 

These  cyclones,  which  closely  resemble  the  Skagerak  cyclones  described  by 
Bjerknes  and  Sclberg  (1922),  generally  form  near  the  southern  end  of  the 
Appalachian  Mountains.  Since  type  C cyclones  were  included  in  the  type  B 
category  by  Miller,  it  was  necessary  to  regroup  the  data  for  1929-1939.  There- 
fore the  frequency  distribution  presented  below  for  this  period  does  not  corres- 
pond exactly  to  Miller's. 

The  complete  tabulation  of  the  three  cyclone  types  for  the  twenty-year 
period  studied  is  given  in  table  1.  The  twenty-year  averages  for  each  month 
and  for  each  year  are  shown  in  figure  45*  A minor  adjustment  for  the  variable 
number  of  days  in  aach  month  was  made  in  computing  the  averages. 

The  variation  of  cyclone  activity  from  year  tc  year  is  shown  clearly  in 
the  table.  In  the  least  active  year  (1946-1947)  only  14  cases  of  cyclogenesis 
occurred  during  the  nine-months  period,  whereas  the  most  active  year  (1948-1949) 
saw  39  cases.  The  average  number  of  cyclones  was  29.9  per  season  of  which  12.6 
were  type  A,  19.8  were  type  B and  5.4  were  type  C. 

It  can  be  seen  from  figure  45  that  September  and  October  are  the  least 
active  months  in  terms  of  cyclone  frequency  with  only  1.8  cyclones  per  month 
compared  with  ,.0  cyclones  per  month  in  February  and  March,  the  most  active 
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Fig.  45.  Monthly  frequencies  of  east  coastal  cyclogenesis. 

Based  on  data  for  the  years  1929  - 1939  and  1942  - 1952. 
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months.  The  seasonal  variation  of  cyclone  frequency  is  due  largely  to  the  vari- 
ation in  the  number  of  type  B cyclones  which  ranges  from  0.3  in  October  to  2.0 
in  March.  The  type  A cyclones  range  in  frequency  from  1.1  per  month  in  October 
to  1.8  per  nont'a  in  November,  while  the  frequencies  of  type  C cyclones  range 
from  0.2  per  month  in  September  to  1.1  per  month  in  April.  More  than  two- thirds 
of  the  B cyclones  occur  in  the  four  months  December  through  March,  whereas  the 
A and  C type  cyclones  shir-  no  such  preference  for  the  winter  months.  It  is  in- 
teresting to  note  that  the  secondary  minimum  of  type  3 cyclone  frequency  in 
February  which  was  discovered  by  Miller,  also  appears  in  the  1942-1952  data. 

A study  was  also  made  of  those  type  B regimes  which  failed  to  produce  type 
B cyclogenesis.  A type  B regime  is  characterized  by  a primary  cyclone  west  of 
the  Appalachian  Mountains,  generally  in  the  Great  Lakes  region,  and  a warm  or 
stationary  front  extending  eastward  across  the  coastal  plain  with  a ridge  of 
high  pressure  to  the  north.  Those  type  B regimes  which  did  not  produce  second- 
ary cyclones  in  the  east  coastal  region  are  here  denoted  type  D situations. 

The  monthly  frequencies  of  type  D cases  are  shown  in  table  2 for  the 
twenty  year  period  investigated.  Also  shown  is  the  ratio  of  number  of  type  B 
cyclones  to  the  sum  of  type  B and  0 cases  for  each  month.  The  last  statistic 
may  be  considered  to  be  a measure  of  the  probability  of  type  B cyclogenesis 
when  the  more  obvious  type  B characteristics  are  present  on  the  surface-  map. 

The  table  shows  clearly  that  not  only  are  September,  October  and  May  months 
of  low  type  B frequency,  but  also  the  probability  of  type  B cyclogenesis  under 
apparently  favorable  conditions  is  low  in  these  months.  Similarly  the  months 
of  high  type  B frequency,  December  and  March,  are  also  the  months  in  which  type 
3 cyclogenesis  is  most  probable  under  these  same  conditions.  On  the  average, 
in  the  nine-months  period  only  thirty  percent  of  the  favorable  situations  pro- 
duce type  B cyclones  in  the  region. 
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TABLE  1.  NUMBER  OF  CASES  OF  CICLOGENESIS 
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